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Exercice l, 

1) Soit y(x) une fonction continue sur [a,b\. Notons x Qt x u ..., x n une suite de points 
distincts de ]a, b\ . 

Notons y(x i% Xj) = 6y[xi,Xj] % y(xo,x u x n ) =é n y[x Q ,x u x n ) 

Montrer que pour tout k 6 {0, 1, ...,n}, on a Va: € [a, b] 

y(x) =y(x< t )-i-{x-xo)y(xo,x 1 ) + ... + {x-xo)(x-x 1 )...{x-x k )y{x t x 0t x u ...,x k ) 

2) On suppose que la fonction y(x) est de classe C" +1 sur [a,b] . 
Montrer que Vi £ [a, b] , 3c x £ [a, b) tel que 

yfr*U (Cj ) 



3) Montrer quey^Cx^o.rr, x„)=y| x,...,x,x ,x u ...,x n | . 

v(p+l) fois 

Exercice 2. Soit / une fonction 4 fois continuement dérivante sur [a, 6] dont on con- 
naît la valeur en des points équidistants xq , x\ , x% , i^ où x^ = a + i/i i = ? 1 , 2, 3. 

1) Donner une valeur approchée à l'aide du polynôme d'interpolation associé à la 
suite (xq,Xi,x 2i Xz) des réels f'(x) où x = x t et x ^ Xi. Expliciter l'erreur commise 
pour chacun des cas. 

2) En déduire que : 

/'(*o) - gj [a/(«o) + bf(xi) + cf(x 2 ) + etf(xs)] + 0ti>fW (U) 

Déterminer a, /?, a, b t c, d. 

3) Soit / : [1,4] -* H telle que f(x) = — considérons la suite SCq = 1*ïi = 2,i2 = 

3,x 3 = 4. X 

3 

Déduire la valeur approchée de /'(l) et /'(-) 

4) Donner une majoration de l'erreur pour chacun des deux cas. 



ETlMiE 



14 



Exercice 3. Démontrer l'estimation du cours pour m = 4, î-fi. si / € C 6 (R, 3t), alors 
il existe une constante C telle que, V h < h^ , 



«M-&2 



/i 4 



<Cfc 8 . 



Exercice 4. Soit / € C 3 (/î, R) donnée, soit x € R et & > donnes. Soit C) = 
a?o + /a, Ta = zq + 2/i et soit 5 la fonction définie par: 

P(x) - /(ro) + ^^(x - *„) + êÊgfA (x - x„)(x - *,)■ 

1) Vérifier que ff(^j) = /(^j) pour j = 0, 1,2 et en déduire qu'il existe £0 € [x ,Xij 
et f, 6 [* x ,«a] tels que /'(&) = /(&) ^ /'(&) = ^(£1). 

2) Soit r la fonction définie par r(x) = /"(x) — 5(1). 

Déduire de 1) qu'il existe rj e [fo, fi ] tel que r"(r]) = et donc r"(i) = / ij*r< s >(t)dt = 

rftIAfc 



/ 



3) Déduire de 2) que [/(x) - gr(»)| < 2/i 3 max l/ (3) (0j si x e fo, ^2] - Comparer 

avec le développement de Taylor, 

Exercice 5. Soit / e C([— 1,-f-l] ,R) donnée et soit p le polynôme de degré 2 qui 
interpole / en les pointe -1, 0,-t-l. Exprimer /*.' p(t)dt en fonction de /(— 1), /(0) 
et /(l). Vérifier que la formule ainsi obtenue coincide avec la formule de Simpson. 
Exercice 6 . 1) Appliquer les règles de Simpson et du trapèze pour n m 6 pour 

calculer l'intégrale de sin x entre et — à partir des 7 valeurs donnés dans le tableau 

suivant: 

x u 12 S 4 3 13 2 

sinx 0.00000 0.25882 0.500O0 0.70711 0.86603 0.96593 1.00000 

2) Comparer les à la valeur exacte.Que peut-on conclure. 
Exercice 7. Le polynôme de Legendre de degré M est défini par 

*•» = WmW & ' l W • 

1) Déterminer £0» £>i /^et £3. 

2) Montrer que pour i,j € N, ï '^ j on a a {L^Lj) = / £j(i)Lj(£) = 0. 

3) Montrer que L M a exactement M zéros réels distincts tous compris dans ]-l, +1[ . 



*-- 



EHMIE 



15 



4) Montrer que la formule de Gauss-Legendre a M points est exacte pour les polynômes 
de degré r = 1M - 1 

Exercice 8. Soit a G R donné tel que < a < 1, soit Zi = -1, t 2 = -a, t 3 = 
Q t c 4 = +1» et soit u>i, W2, ws, W4 quatre nombres réels. Nous considérons la formule 

4 

de quadraure définie par J(g) =£ ^g(tj) t où g Ç C([-l,+l]) donnée. 

j=\ 

/-t-i 
p((t)dt, pour tout 

polynôme p de degré '6. 

2) Existe-t-il a tel que J(p) = / p{t)dt, pour tout polynôme p de degré r, avec 

r>3? - 7 - 1 

Si oui, quelle est la valeur maximale de r et que valent aloxs a et Ui , u^, oj 3 , uj 4 ? 

Exercice 9. Soient t, = -1, ^ = 0, t 3 = a, où a € ]0, 1] , trois points distincts fixés 

dans rintervaUe [-1,1]. 

Si g est une fonction continue définie sur [-1, -+-1) , on considère la formule de quadra- 

d 

ture à trois points J(g) ^S^ ^M^)- 

où ijjj, j= 1, 2, 3 sont les poids de la formule de quadrature. 

1) Calculer les poids us it u^ t u/ 3 en fonction de a pour que la formule de quadrature 
soit exacte pour des polynômes de degré 2. 

2) Existe-t-il a pour que la formule de quadrature soit exacte pour des polynômes de 

degré 3? Si oui, donner la valeur de a et les poids u x , u 2 , a; 3 . Sinon, justifier votre 
réponse. 

3) Existe-t-il a pour que la formule de quadrature soit exacte pour des polynômes de 
degré 4? Si oui, donner la valeur de a et les poides wi, a^, u/ 3 . Sinon, justifier votre 
réponse. 

Exercice J.Q, Soit -1 < t t < t 2 < +1. Etant donné ces deux points, nous cherchons 
deux poids w u o/ a qui définiront la formule de quadrature: J(g) = wtffo) 4 u) 2 g{t 2 ), 

pour approcher / g(t)dt, la fonction g étant continue sur [-1, +1] 

-î 

1) Trouver les poids u u un en fonction de îi, t 2 tels que J(p) = f p(t)dt pour tout 
polynôme p de degré 1. 

2) Que deviennent les poids lorsque t\ = -f 2 

3)On considère le cas où ^ = -1 et t 2 = +1- La formule est-elle exacte pour des 
polynômes de degré 2 
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4) On considère le cas où t 3 - -et et t 2 = +a avec < c* < 1. Pour quelle valeur de 
a, la formule de quadrature est-elle exacte pour des polynômes de degré 2? 
la formule de quadrature ainsi définie est-elle exacte pour des polynômes de degré 3? 
la formule de quadrature ainsi défini est-elle exacte pour des polynômes de degré 4? 

Exercice 11. On définit la fonction g(x) = Var € J-co, 1[ 



1) Calculer F(x) = f g[t)dt 



X<1 




2 

2) Quelle est la valeur de F(x) en x = - 

3) Donner le degré et l'expression du polynôme de Lagrange qui interpole g(x) aux 
points 0,-,-. 

4) Trouver les coefficients 00,01,02 tels que pour tout polynôme P de degré < 2 on 
ait: 



I P{x) dx = coP(O) + ciP(l) + c 2 P(2) 



5) En utilisant un changement de -variable, déduire de la formule précédente les coef- 
ficients do,d\,d2 tel que pour tout polynôme Q de degré < 2 on ait: 



a 

ï 

/ 



Q(x) dx = d«Q (0) + d x Q (j\ + daQ (| 



Utiliser cette formule pour donner une valeur approchée de log 3 

Exercice 12. Soita e ]0, 1] C R donné, soit t x = -a, fe = 0, £ 3 = a, et soit 

uj\, a»2, wj, trois nombres réels. Nous considérons la formule de quadrature définie 

par J{g) =£ ujj g(tj), 

où jeC{[-l t -fl];R) donnée. 

p{t)dt, pour tout 
polynôme p de degré 2. 

p(t)dt pour tout polynôme p de degré 
■6. 

3) Existe- t-il a tel que la formule de quadrature soit exacte pour les polynômes de 
degré 5? Si oui, calculer a et comparer avec les zéros du polynôme de Legendre de 
degré 3 
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Série n°2 

Correction d es Exercices 
Exercicel: 

y{x) est une fonction continue aur [a, b] , montrons que pour tout k € {0, 1, n} on 

CtaiT VM * {X ~ xM *°' Xl) + - + {x " Xo){x ~ *>>-<« - qOvfe**, ...,**) 

y( x ) = V(xQ) + (x-x Q )y( Xi x ) 

y(x,x ) = tf ( Zo , xO + (x - a^yfoafc,**) 

»( I .a;o,*i)=îr(xo,a:i,aî2)-l-{a;-ara)3/(r,xo l * ll x a ) 



V(*.*<>.fi .-•^-i) = !/(xo 1 x Ij x 2 »t) + (*-x*)v(«,«tt,*i,ï2 X k ) 

On multiplie la seconde équation par(x - x ) la troisième équation par (x - x ) (x - Xi) 
et ainsi de suite la dernière équation étant multipliée par (x - mù(9 - x, ) (x - 

x*-i),on obtient: ' v 

y( x )=y{xo) + (x-x )y(x t x<i) 



(x - *o)(x - x0 (i - efc-Ovfeao,^ x^, « fc _ %)(â] _ 3 , . _ 

nn^ Û,X,, ^ , *;' ; f t)+(x " Xo,(x "" Xl) (*-s*-i)(s-z0î/(x,* ,s,,t 2 *) 

On additionne le tout et après simplification on obtient: 

»(*) = y(^o)+(x-xo)y(xo,xj)+ + (*-*q)(x-xi).... (x-xOii(x r« r x ! 

2) On a d'après la question précédente: * ** ( ' °' b "'^ 
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y(x) = P n (x) 4- 7r(x)y (x, x , x, , , x n ) où i> n (x) est le polynôme d'interpolation de 

Newton associé aux points distincts xo.xi, x n etn(x) m (x - x ) (x - Xi ) (r ar, 

Or on sait que Vx € [a, 6] , 3 c* € [a, 6] tel que : 

y(x) = P n (x) + tt (x) T^y,^ 11 <*) 



-) 



rv s ( \ V ( " +1) (*) 

Doù : y(x,x ,Xi,...,x„) = . ^*n 

3) Montrons par récurrenceque: 

yW(x,xo f Xi, ,x n ) = y I x.J.^.^ Xo^i, X 

d u(g + /i,x ,si,...,= Cf») -vtx.xo^x^.-.xj 
-y(s,x ,x : ,...,x n ) =bm £ 

Puisque la différence relative est symétrique 2 à 2 on obtient: 
d , , .. y(x + /t 1 x ,x 1 ,...,x n )-y(xo,xi,...,Xmg) 

-y(x l x ,x 1 ,...,x„)=hrn g 

— y(i,x ,x lf .. M x„) =hmy(x-\-h i xo,Xu.. n x n ,x) 

AI n—0 

— -y(x,x |Xii-.Xn) =lim y(x+ fc.x,x »*ii- »«n) 

—•y(x,iD.Xi,... t x n ) = y{tft<ttfAtof*U'«»fti) 

Supposons que la propriété est vTaie jusqu'à Tordre pcà.d y (j?) (x.xoii'i, ,x„) 

y (x,! e,«o,aîi, x„ I et montrons la à l'ordre p+1 



*v 



ytP+1) (a;, xo, Xl , x n ) = £y (p > (s, xo, H _ In) 

^^(x.XO.Il, .^rt) = T'y I g^.""-g |gO)gli »n 

\ (p+lî/oH 

^^(x.xcxj ,i„) = y I x,x x,x ,x u x n J d'où le résultat. 

Exercice 2: 

1) Cherchons le polynôme d'interpolation de Lagrange P^ (x) de degré 3 associé â 

/ aux points Xo, afe, 2:3, 13 où x t = a + «A, i = 0, 1, 2, 3 

Soit, (yo, <pi,<p2) <Ps) la base de Lagrange associée à ces points. 

On a P 3 (x) = £>< (x) /(x,) où 
1=0 
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Vi(x) = 2/i 5 < x ~ x °) (* " l2 ^ ( x ~ X3 ^ 
Vato = -53 ( x - x °) ( x " x i) ( x " x ») 

VsW ■ 53 C*'~ x °) ( œ " Xl ^ < x " * a ) 

Or on sait que: /(a;) - P 3 (x) + E{x) où £(x) - -^ (x)/ (4) tf«)i avec »(#) = 

.=3 

Jl (s -x0 et( x Ç [x 0l x a ] 

donc /'(x) = r^(x) + £'(x), d'où la valeur approchée du réel /'(x) c'est P&x) avec 
une erreur de E' (x) 

D'où: si x = x„ /'(x,) - 2 4 W /W + ^ W /*•&*) 

si x 9fc x„ /'(x) = gv>S (x) /(x,) + ^ (x) /W(fc) + ^x (x) /»&} 

2) si x = x f( XQ ) = jg^j M /(**) + J^ W / (4, «-o) 

Cherchons <^(x ) pour ; « 0, 1, 2, 3 et t* (x ) 

<4>M = " gj-g K* - *») (* - «*) + (* - **) (* " x a) + (x - xi) (x ~ s»)] 

#(*) = ^3 [(x - x 2 ) (* - ia) + (* - *û) (a ~ »»3 + (■ " x o) (x - *a)] 
^i(x) = --^ [(x - xt) (x - x 3 ) + (x - xo) (x - x a ) + {x - x ) (x - a*)] 

v 3 (x) = ê^â K* - «aJ (* ~ x =) + ( r - x °) ( x " Xl ) + (:r " Xo) (x ~ Ia) l 

D'où ^ (xo) = ~, ^i(xo) - |, ^(x ) = -^, tfi(ao) - gr, *' (*>) = -6/l 9 
Donc /'(xo) = ^ [-ll/(xo) + 18/(x a ) - 9/(12) + 2/(x 3 )j - -tf/W fcj 

D'où le résultat avec a - 6, a = -11, 6 = 18, c = -9, d s 2 et = -- 

Exercice 3. / € C 4+a (R, R) = C^R, R) (m = 4). 
D'après le cours , nous avons 



et 



€(/W) - /ta + a) - VW + /ta - *>, 
Biffa) = fc(tf M) - WfcfflW"*» - s Mf (*<>))■ 
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Puisque l'opérateur fijj est linéaire et en utilisant l'égalité ci-dessus, nous obtenons 
6JJ/(*o) = Slf{x -h h) - 2#/(*o) + èlf(x - h) 

( or par définition 8 h f{x ) = f(x + j) - /(ar - -)) donc nous avons 

Biffai) = /(io + 2/i)-2/(io-h/i)+/(ïo)-2/(xo + ^) + 4 /( :E o) 
-2/(i - h) + /(io) - 2/(xo - h) + /(r - 2k) 



D'autre part, le D.L à l'ordre 6 de / au voisinage de i nous assure que 



V 



/(xo-2*) - /(x )-/'(xo)2n + r( I o) i ^-/ f3) (xo)^ L 



+ 



+rw ç_ /(sl(I0) ç +/Mw ^) 



6! ' 
/(*. + h) = /(*<,) + f {*>)h + /"(xo)| + / <3) (^o)| + / l4) (*o)^ + /< 5> (*o)§ + 

/(xo - fc) = /(«.) - /W + f(«o)| - /«WÏ + / (4, (^o)^ - /< 5) (*o)£ + 

où (iJi,»h 1 »fi,'?4)G ]sso,*a+2&[ x]x -2/i,ro[x lx ,i-o + t| * }x - h,x [. 
Après substitution dans la formule de 6£/(xo), nous obtenons (6J/(io) = /(io+2A) + 
/(xo - 2fc) - 4/(x + fc) - 4/(x - h) + 6/(x )) 

«/(*■>) - / (4) (x )fc' + fê(/<%i) + / (8, (%)) - ^(/ (6, (*.) + / w (*)>] Z» 6 - 



Soit hc € ]0, +oo[ fixé et soit 
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C = — r max 

90 i -2ho<t<*+'iAio 



|/ (6, W| 
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Exercice4: 

/ e C 3 (H, fi) donnée , soit x «G fi et h > 0, Xi = Xo + h, xa = xq 4- /i 

*(*) = /<*o) + ^i(* - *o) + %^(* - x )(x - wù 

où A h /(x ) = f{x -h h)- /(x ) et / e C 3 (R, R), x c = x + i/i (i = 1, 2). 

1) Il est clair que g(xo) = /(xo). D'autre part, nous avons g(xi) = /(xo) + Aj,/(xq) ™ 
f{x + ft) = /(xi), et 

£?(* 2 ) = /(xo) + 2A ft /(x ) + A=/(x ) 
9fa) = 2/(x + A) » /(xo) + A2/(x„). 

A?/(x ) = A,,(/(x a + *) - /(x )) = A„/(x + h) - A„/(x ) 
ûï/(*o) = /(xo + 2/i) - /(x & + A) - /(x + ft) + /(x ) 
^/(*o) = /(x + 2&) - 2/(x -h /*) + /(x ), 

donc p(x 2 ) = /(xo + 2ft) = /(x 2 ). Posons r(x) = /(x) -y(x). 

Puisque r(x ) = r( Xl ) = r(x 2 ) = et puisque r e C^K.K), on peut utiliser le 

théorème de Rolle pour obtenir 

3Cq e ]x ût xi [ tel que r'(Ço) = i.e. /'(Ço) = g'(Ç Q ), 
3Ci e I*I,a*| tel que r*(Ç,) = i.e. /'(£) = /(Ci), 

2) Puisque r'(Ç ) = r'(&) ■ et puisque r e C^fR,!*) on peut à nouveau utiliser le 
théorème de Rolle pour obtenir 3rç e ]Ço, Ci[ tel que r"(n) = 0. 

Par conséquent, puisque r'" est continue, nous avons r"(x) = r"(x)-r"(r}) = / ;f r< 3 )(0<fc 

Puisque £ est un polynôme de degré deux il est clair que r< 3 )(f) = /®(t) - 5 W (') = 

/< 3 >(0 et donc r"(x) = [ rf f™{t)dt. 

3) Considérons par exemple le cas où x e [x ,xij . 
Le cas où x € [x^, xj\ se traite de manière analogue. 

D'après 1) , r(x c ) = et donc /(x) - g(x) = r(x) = r{x) - r(x ) = f r'{s)d3. 

Par conséquent, |/(x) - y(x)| < / |r'(s)|rfs< max |r'(s) ; . (+) 

soit s 6 [x ,Xi] - D'après (a), r'{Ç ) = et donc r'(s) = r'(a) - r'(Ço) = / r"{t)dt 



Par conséquent, \r'(s)\ < f \r"{t)\dt\ < f ' |/(i)| dt < h m 



max r"(c)L («* 
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Soit t € [a;o,sri] ■ D'paèrs (b) r"(t) = f rr7 (3ï (") du > & P^ conséquent 



r"(t)\ < 



J3 

f TffW(u)du < f |/ (3) W| du < 2h max |r (3) (u)| . (***) 

«g 



Les inégalités (*), (**) et (* * *) impliquent donc 

I/O») - *(*)l < 2A 3 Jg^ ir< SÏ (ii)| ■ 

Comparons ce résultat avec la formule de Ttoylor. 

Soit i e [îfo.Xî] donné, puisque / € <7 3 (R,R),3Ç 6 [ar 0ï s] tel que 

/(x) = G(x) + ^^tr-*o) 3 
où G est le polynôme de degré deux (défini par 

G(x) = /(xo) + /'(*«)(* - x ) + ^y^(* - *o) 2 - 
Par conséquent, lorsque x € [x - x Q ] , nous avons 



|/(x)-G(x)| = 



/ (3 >(0 



[x - x ) 



< &-**> max |/W(t)|, 



Interprétation; Nous avons donc à notre disposition deux polynômes de degré deux, 
g et G , permettant d'approcher une fonction / € C 3 (R,R), au voisinage d'un point 
xq S [x 0t x 2 + 2h\ Terreur maximale entre / et les polynômes g et G est d'ordre trois 
en h sur [x ,x + 2h] . Le polynôme G fait appel aux dérivées première et seconde de 
/ au point x . Par contre, le polynôme y fait appel à des approximations numériques 
de ces dérivées. 

Exercice5: 

Soit / e G ([-1,+1] ,R) , alors le polynôme P de degré 2 qui interpole / en les points 
t = — l r *i — et t 2 = +1 s'écrit ainsi 

. P(«) = /(-!)«.(*) + /(O)Vi(t) + fWf2(t) 

où <pu, y>i, ip2 est la base de Lagrange de ïï> 2 associées aux points to = —l» *i — 
0, ^ = 1 et est explicitée dans le cours exemple 1.1 (chl). i.e. 
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un calcul simple donne 

+1 +1 +1 +i 

y , j<t)(ft=/(-i)y"w(*)A+/(o)y^(o*+/(i)/ww* 






/-m y f*' 4 [■ ■ 1 

^o(É)dfc=-, / ^îW^t = -, / <P2{t)dt — -,et par conséquent 
-i 3 J -i 3 J -i 3 



-î 



/pW««=|(/(-i)+i/(0) + /(i)). 

-1 
Il semble donc naturel d'approcher / f{t)dt par la quantité */(/) définie par 



/(/) =-§(/(-!) + 4/(0) + /(!)). 



Dans le chapitre 3, on appelera y(/) formule de quadrature. Par construction, </(/) 

*+l 

intègre exactement les polynQines de degré 2 au sens où / q{t)dt = J(q) Vç 6 OV 
Cette fonnulc de quadrature s'appellera formule de Simpson. 



Exercice6: 

On va appliquer les règles de Simpson et du trapèze pour n = 6 pour calculer l'inté- 
grale de ain x entre et — à partir des sept valeurs tabulées en la table suivante: 

X " 12 6 4 I 12 2 

sini 0.00000 0.25882 0.50000 0.70711 0.86603 0.96593 1.00000 

Après calcul la règle du lrapèze a donné 0, 99429 et celle de Simpson 1, 00003 
Eu comparant à la valeur exacte soit 1 on obtient: 

|1 -/ri =3 5,71. HT 3 

|1 - /ttmpl * 3 10- 5 

On voit très bien que la règle de Simpson est meilleure. 

■ 

Exercice 7. Soit M g H, L M (t) = ^^ [(*» - 1)*) 
1)M0 = 1, W) - ~((* - V) = t. 
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(f*-i)*+A<*-i) a 






4x2 



i 



MO = I {3t 2 - 1 + 2t 2 - 2} = ^(5< 2 - 3) 
MO - 1'(' 2 - 1) 

2) Montrons que (X», L,) = I L x (t)Lj(t) = si t £ j, en effet, supposons 
obtient alors en intégrant par partie 



i > j; on 



-i 



-i 



Puisque (< 2 - 1)* a un zéro d'ordre i en 1 et en -1, la (i - l) ,èmc dérivée de {t 2 - 1)' 
s'annule en / = 1 et en t = — 1. Ainsi nous obtenons 



-i 



-i 



En intégrant par partie j fois comme ci-dessus, nous obtenons: 
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3) Soit t u t 2 , ...,t a les points strictement compris entre -1 et + 1 en lesquels L M 
change de signe. 

Clairement ces points seront des zéros de L M et on a donc s < M. Si on pose 

p(t) = {t~ ti)(t - i 2 ) (f - t s ) alors p e P, et, puisque p change aussi de signe en les 

points U % l<j<soa obtient p(t)L M {t) > Vt 6 [-1. + 1] ou p(t)L M (t) < V* € 

i-l ) +li.Danstouslescas î puîsquep(<)/ (M (i)^0 1 olia (p.Ijf) ^ 0. 

Puisque L Q J n ,.„,L M forment une base de P M alors il existe a ,a It ,», tels que 

P(0 "/) %%(*) et en utilisant ia partie 2) on obtient 



1=0 






(p, L M ) =J2 <** / h{t)L M {t)dt = a, (L a ,L w ) . 
fi* i x 

p{t)L M (t)dt ï on a nécessairement s = M et donc les M 
zéros de L M sont ^ f fa, ...,* A/ . 
4) Nous aUons montrer que la formule de Gauss-Legendre a M points est exacte pour 

les polynômes de degré r = 1M - i. Pour cela, considérons J(g) =JT ^g^) la 

formule de Gauss-Legendre à M points et soit p un polynôme de degré^À/ - L alors 
nous pouvons définir pour t € R 

où fcu**,...,*? W est !a base de Lagrange de P M _ : associée aux points de Gauss 
Autrement dit, p est donc l'interpolant de p de degré M - 1 aux points de Gauss 

Soitç(f)=p(<)-p(i) Vi€R,aJors rfe 5 (ç) = 2M-1 tels que q(t : ) = 0, j = l M 
Ainsi g est divisible par le polynôme v de degré M défini par; u(£) = (f - /,)(* _ 
ï 2 ) ; ....(i - t M ) WeR c-à-d il existe un polynôme w de degré Af - 1 tel Que 
</(*; = v(0u>(ï:) Vf € R. 

Puisque v est un polynôme de degré M qui s'annule en les M zéros de L M qui est lui 
même un polynôme de degré M, 3a e R tel que v(t) = *£*(*) Vt € R. 
Puisque w est un polynôme do degré M - 1, il existe &, A, ..., M ^ tels que 

M-l 

w(t) «£ &!,(*). 

r '0 
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Ainsi donc,en utilisant la propriété d'orthogonalité des L ? (t) (voir 2)),nous avons 

I q(t)dt = / v(t)w(t)dt ma £ & / Lu(t)L k (t) = 0, 
-i -i H> I, 

donc nous avons prouvé que 

-fi +i 



J Vif)dt ^ j p(t)dt, 



-î -i 

et en fin, par définition de p, nous obtenons 

+î U +1 






ce qui achève la démonstration. 

Exercice 8. 

1) D'après la formule du théorème 3.2 du cours, nous obtenons 



+i 



où les fonctions tp h j = 1, 2, 3 et 4 sont les polynômes de la base do Ugrange de P, 
associée aux points U = -1, t 2 = -a, t 3 = a ei U = +1. 

Puisque çtfft «fi ff "7 (Ici n = 3), alors nous avons 

* + Q *-tt t - 1 .. t -f 1 £ - l « - i 



/,» viu * — u t — 1 



-a -f 1 — a — a ' -o — 1 

et donc 

+i 



+i 



_I -1 

Pour des raisons de symétrie, u A = Wl et^=^. Par construction J(p) intègre 
exactement tout polynôme de degré 3. 
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2) Tout polynôme de degré r peut s'écrire p(t) = at" + 9 («) où g est un polynôme de 
degré r - 1 et a e R. Par conséquent, 

4 

7(p) = a ]T w$ + J( g ) et 



+i +i +i 

/ p(t)dt = a ft r dt + fq(t)dt. 
-i -î _i 

Donc, pour que J[p) = H p(t)dt pour tout plynôme p de degré r, il suffit que 

'^ = 7 -, *^* P0Ur t0Ut P 01 ^ ™ 9 de degré r - 1 et que J(&) = (^ fdt. 

Nous procédons donc par étapes pour déterminer le degré maximal du polynôme pour 
lequel la formule de quadrature est exacte. 

wZ£V / QUe ^ f0nnUle ^ QUadratUre GSt «K"* 5 P our les Polynômes de degré 

^ t*dt on obtient o = _ et ainsi pour cette valeur de a, la formule de 
quadrature est exacte pour tout polynôme de degré 4 et nous obtenons «* = W « •= I 
et wa = u 3 = g et la formule de quadrature s'écrit: J(g) = - \g{-l) + g{l)\ + 

*J£Z!? e8t faCilC ^ VérLfîer QUe œtte f0rmUle ** aussi ***** P our te Polynôme 

0= ft*dt 

ïlle ne l'est, plus pour le ni vu 
En effet, 



t puisque 

Par contre, elle ne l'est plus pour le plynôme t 6 . 



^)=i((-D 6 + (i)vi((-^)^(^)=i 
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Ainsi la réponse à la 2 timc question de l'exercice est positive et on a r — 5. 

Exercice 9: 

1) D'après la formule du théorème 3.2 du cours, nous obtenons 

+ i 
w 3 = J nWdt, 

où les fonctions <fj, j = l,2,et 3 sont, les polynômes de la base de Lagrange de Pa 
associée aux points Éj = — l, 1% ~ et £3 = a 

Puisque <fj(t) = n 77 7^ (Ici n = 2), alors nous avons 



tt - t 2 )(t - t 3i 



t(t - a) 



n M t-ttXt-t») (t + l)(t - a) 

„ Ml = ( ( -«!)(«-«») = (' + 1 ) f 
m> (*l-tl)(«I-tï) "(« + !)« 

/•+! 1 /•+> 1 



3 



o£^ 

2 



+1 



i 



l + o 3 



w 2 



3(1 + a) 
.i£V + (l-a)*-a— I[| + (i-a)f-«t 

1 fi l-a 



+1 



J -1 



**•— â 3 



2 

+1 



1 (1-tt) 

ar + j ^ « 



TD 3 = 



o(q + 1 / -1 



o(o+l) 



f 2 +f = 



1 



û(a+l) [3 2 



1-W 



, -1 
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3 + 2 + 3 2 



3a(a+l) 



2)Par construction, /(p) intègre exactement tout polynôme p de degré 2. 

r-H 

Pour que J{p) = J p(t)dt pour tout p G P 3l il suffit de vérifier que J{t*) = 
/ t 3 dt; en effet 

J(t 3 ) = -W! + wga 3 
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\ 



/+1 
Pât = j [t^Vx = 0. 

donc Jfr 3 ) = -wi + W3Û 3 = 

Q 2 = let a^-l=»a= 1 

donc wj = iu 3 = - et ït*2 = « d'où on obtient la formule de quadrature de Simpson. 
3 . 3 

/1 2 2 

* 4 d* = r -/(O = ^ donc il n'exixte pas de a Pour que la 
_i o J 

formule de quadrature soit exacte pour les polynômes de degré 4. 

Exercice 10; 

Soit -1 <rj <t 2 <let J{g)^uj l giU)+uJ29{t2) 

1) Cherrions les poids bft,<ife en fonction de £1 et t% tel que J{p) — I p(t) dt pour 

-1 

tout polynôme p de degré 1 Les valeurs de oJj et W2 en fonction de t\ et £2 tel que 

J(p) =^2 ^jp(^) = / p(*)<if pour tout pf?i sont données par 



'i = I m^ 



■1 

où les fonctions <^,, i = 1,2 sont les polynômes de la base de Lagrange de Pi associée 

aux points <j, U Puisque <f*(t) = U 7 y (ici n = 1). Alors nous avons 

■-1 (^ - fi) 



■I e ; 






"H 



ipi(t)dt = —, 

wa = / «^(O*"- - : f- 

J -1 *2 - n 

2) Si ti = —£2 alors u;i = u>2 = 1. 

3) Si *i = -1 et Î2 = 1 alors u»i = \a% = 1, 

/ +1 I 2 

**<** - 3 I< 3 ]-î = 3? or */{t a ) = 2 donc la formée n'est pas exacte 

pour les polynôme de degré 2 

4) Si t] = — o et ti = a alors Uj = u;a = 1 
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cherchons a pour que la formule de quadrature soit exacte pour les polynômes de 
degré2 J(fi) = / Aft -+ («, +WS )rf = | Comme < o < 1. on a J(«) = 

J _P{t)& Vp G 1%, si et seulement si a = -L La formule de quadrature Jfo) « 

* 9 l~7!/ + 5 [ vf / ^ d ° nC exaCte p0Ur IeS P ol y nôme d « «Nrt r = 2. De plus 
, lorsque p(e) = f», n a J(p) - ^ p ( t )dt = 0, la formule de quadrature J(g) ~ 

9 \"y/3j + 9 X'js) CSt GXaCte pOUr deS P o| y nôme de degré r = 3. On vérifie sans 
peine^que cela n'est pas le cas pour des degrés r plus élevés en prenant par exemple 

Exercice 11: 

On définit la fonction g(x) = — î— V* G )-oo, 1[ 
l)F(x» = -log(l-*) l ~~* 

2) F(|) = log3 

3) Le polynôme de lagrange qui interpole g en trois points 0, |, | est de degré 2 et il 
est donné par: P a {x) = |i a + 1 

4) On pose p(x) = l, p (a:) = * et p(x) = a» on trouve le système suivant: 

y ldr«2=Cb + Ci +C2 



o 
■» 



/ 



Jx*dx~ g-^4-4fe 



x dx = 2 = c, + 2ca -=> co = cj = i Cl = 1 
- 3 3 



5) On peut effectuer le changement de variable x = 3y on a: / p{x) dx . 3 7 ? ( y ) <ty 

1 , - o o 

ce qui dorme <*, = -q d'où d = d 2 = -, i/i = - 

Utilisons cette formule pour donner m?e valeur%prochee de log3 
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• 


2 

2 \ 

^3) = log3 = y g(y)dy 

2 '2 

3 3 

Or 1 g{y)dyesi approchée par / 


• 


* 




on trouve Iog3 fe — . 




Exercice 12 



l)Soit a e ]0, 1] c R. Les valeurs de u; t en fonction a tel que /(p) =V] w >0&) ■ 

r 1 

/ p(*)di pour tout p 6 1*2 sont données par 






41 

-1 

où les fonctions gfo t = 1, 2, 3 sont les polynômes de la base de Lagrange de P a associée 
aux points fc = -« : i 2 = 0, t z = tt . Puisque ^(t) » n ~^r (Ici n = 2). Alors 

nous avons 

/,* i i-a i+at-a 

^ l() "^'^T' ^' = • 

o zq a — a 



donc 

2a 2 ' 






/ +1 -2 



Pour des raisons de symétrie, c* ~ «*> Par construction, J(p) intègre exactement 
tout polynôme p de degré 2. 

2) Pour que J{p) = J p {t)dt pour tout p € J> 3| il suffit de vérifier que J(t*) = 
/ r 3 ^; en effet 

J{t 3 ) = wi (-a) 3 + tJa.O + u; 3 a 3 = 0, 
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+1 i r+i 



et, d'autre part f t*dt = i [t 4 ]+\ f Pdt = ~ [t 4 ]!î = 0. 

3) Nous avons J(p) = I p(t)dt pour tout p 6 P 4; s» ^(* 4 ) = / l * dt - Puisque 

f+ l 2 

J(f 4 ) = wi(— o) 4 + u^.O + w 3 a 4 = 2a>itt\ et puisque / ( 4 df = -, il suffit donc que 

Wi q 4 = - c'est à dire a = J- e ]0, 1] . 

U V M 

Finalement nous vérifions facilement que 7(t 5 ) = / t 5 dt = et donc la formule do 

quadrature est exacte pour un polynôme de degré 5 lorsque a = J-. 

5 3 

D'après l'exercice 7. 1/, le polynôme de Legendre L3 est donne par L$(t) — -zt{t 2 — -), 

2 

il s'annulle donc pout t = et t - ±\j~* ^a formule de quadrature correspondant 

au choix a = J - est donc la formule de Gauss-Legendre à 3 points. Elle est bien 
exacte pour les polynôme de degré r = 2 ^ 3 — 1 =5. 
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